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Abstract 



The North Atlantic Oscillation (NAO) monthly index is studied from 1825 
till 2002 in order to identify the scaling ranges of its fluctuations upon dif- 
ferent delay times and to find out whether or not it can be regarded as a 
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ation analysis both indicate the existence of weakly persistent long range 
time correlations for the whole scaling range and time span hereby stud- 
ied. Such correlations are similar to Brownian fluctuations. The Fokker- 
Planck equation is derived and Kramers-Moyal coefficients estimated from 
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the data. They are interpreted in terms of a drift and a diffusion coeffi- 
cient as in fluid mechanics. All partial distribution functions of the NAO 
monthly index fluctuations have a form close to a Gaussian, for all time 
lags, in agreement with the findings of the scaling analyses. This indi- 
cates the lack of predictive power of the present NAO monthly index. Yet 
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there are some deviations for large (and thus rare) events. Whence sugges- 
tions for other measurements are made if some improved predictability of 
the weather /climate in the North Atlantic is of interest. The subsequent 
Langevin equation of the NAO signal fluctuations is explicitly written in 
terms of the diffusion and drift parameters, and a characteristic time scale 
for these is given in appendix. 

PACS numbers: 05.45.Tp, 05.45.Gg, 93.30.Fd, 89.69.+x; 02.50.Le, 05.40.-a, 47.27.Ak, 
37.23. Ge 



I. INTRODUCTION 

In order to establish a sound understanding for any scientific phenomenon, one has 
to record numerical data and from the latter to obtain laws which can be next derived 
theoretically from so called first principle models. The so called inverse model method, 
starting from raw data and using statistical analysis as a first step, is of great interest 
since it is model free. Some difficulty arises in particular in nonlinear dynamical systems 
because of the need to sort out noise from both chaos and deterministic components 1,2 . 
Whence to extract meaningful model-free dynamical equations from chaotic-like data is 
an enormous challenge 3 . Practically one is often led to empirical relationships. This is 
often the case in the meteorology/climatology field where there is a widely mixed set of 
various (sometimes) unknown influences, over different time and space scales. Often the 
fast variations are taken as noise terms in some sort of Langevin equation(s) 4,5 as for the 
el Nino Southern Oscillation Index (SOI). 

In order to quantify weather and climate events in Europe and report large-scale vari- 
ability an index has been imagined 6 the so called North Atlantic Oscillation (NAO) in- 
dex (http : / /www.ldeo. columbia.edu/NAO/; http : / /www. met. rdg.ac.uk/cag/ NAO/; 
http : / / www.cru.uea. ac.uk/cru/info/nao/; 

http : //www — bprc.mps.ohio — state.edu/gpl/NAO/Naobibliography.htm). It is 
the normalized sea level pressure (SLP) difference between a station at Ponta Delgada, 
Azores and one at Akureyri, Iceland. 

Since about the mid-50 's the NAO index was trending from negative to positive values, 
but is mostly positive since 1980, a variation attributed to global warming. It is thought 7 
that the influence of slow changes in the ocean and in the greenhouse gases maybe picked 



up as the fundamental causes of a prolonged (upward) trend. 

Until a few years ago, the NAO was not receiving intense attention 7 , because it was 
thought that its phase and amplitude were rather unpredictable, because both involve 
many (time and space) scales which are often intrinsic to chaotic behavior; see also 
reviews 8 10 . Yet, evidence has been presented that NAO exhibits 'long-range' dependence 
having winter values residually correlated over many years, with short-term 2-5 year 
variations and decadal trends. Note that Wallace 11 has argued that the NAO is a local 
expression for a Northern Annular Mode (NAM), also called Arctic Oscillation 12,13 . 

In view of the above it seems of pertinent interest to consider again the NAO and 
adopt specific data analysis techniques when searching for scaling ranges and stochasticity 
features. We start with the Hurst (R/S) method 14 ' 15 followed by a detrended fluctuation 
analysis (DFA) 16 ' 17 of the monthly averaged NAO signal. Such tests supplement classical 
analyses based on frequency spectra 18-20 which are debatable due to the non stationarity 
of the data. Interestingly the data histogram have so called fat tails, resembling the Levy 
flights, signatures of self — organizing systems, today emerging in many areas of physics 
as those mentioned here above. Again these facts seem to exclude low dimensional chaos 
but support the conjecture of Markov dynamics for atmospheric evolution, as already 
suggested in fact many years ago 8 ' 21-23 . 

In the following, we adopt a Markov assumption in order to derive the FPE and 
to write down the Chapman-Kolmogorov equation for the conditional probability of the 
increments Ax (of the NAO index) over different time intervals At. This leads to a 
numerical derivation of the Kramers-Moyal coefficients which are the moments of such 
probability distributions. Up to the second moments, this leads to the diffusion and 
drift coefficients appearing in the Fokker-Planck equation (FPE) and are basic to the 



Langevin equation 24 . It will be noticed that the analytical form of both drift D^ and 
diffusion D^ coefficients are simple. It will be found that the experimental probability 
density functions (pdf ) have all a Gaussian form when excluding the (rare) large (so called 
extreme) events. 

The methods, applied in this paper, are briefly explained in Sect. 2.2 and 2.3 : (1) 
the rescaled range analysis and (2) the detrended fluctuation analysis. In all cases results 
are tested against surrogate or shortened data for error bar evaluation. Both methods 
lead to an exponent characteristic of the classical random walker position fluctuation 
correlations. Next, in Sect. 3, we examine how to describe the statistical evolution of 
increments for different time scales i.e. establishing a Fokker-Planck equation within a 
Markov process assumption. A few comments pertain to considerations on NA weather 
causes and predictions in a discussion and conclusion sections. 

II. DATA AND THEORETICAL ANALYSIS 

A. Data 

The 
monthly averaged NAO index (available on the web sites http : / /www. cru.uea.ac.uk/ 
ftpdata/nao.dat or http : // www.cru.uea.ac.uk/cru/data/nao.htm and updated at 
http : // www.cru.uea.ac.uk/timo/projpages/nao_update. htm), i.e. as the normal- 
ized sea level pressure (SLP) difference between a station at Ponta Delgada, Azores 
(26°W,38°N) and one at Akureyri, Iceland (18°W,66°N) is represented on Fig.l from 



January 1825 to November 2002 (2135 points). 1 

It is a standard procedure that in order to reduce spurious noise effects, the study is 
performed on the integrated series (Fig. 2). Such values can be interpreted as mimicking 
the successive positions of a random walker 15 . The amplitude correlations should allow 
us below to understand the drift and diffusion process (as that of a walker). 

B. The rescaled range analysis 

Introduced by Hurst 14,29 , the rescaled range analysis method computes a ratio R/S 
defined as follows. The time series X = {x t ,t = 1,...,N} is divided into / intervals 
of equal length n. In the k th box, (k = 1, ...,/), there are n elements, XJ (n) = 
{xj,j = (k — l)n + 1, ..., (k — l)n + n{= kn)}. The local fluctuation at point j in the k th - 
box, i.e. (xj — {xy n ') is calculated as the deviation from the mean {xy n ' = ^ YJj=i x \ , 
in that A; t/l -box. The cumulative departure Y^\n) up to the m th point in the k th -hox. 
(of size n) is next calculated 

m m 

Yi k \n) = £ (xf - (x)i fc) ) = (£ 4 ] ) ~ m (x)^ (1) 

for m — 1, ...,n and in all k boxes and where (x)„ = - I]™_i x) . The rescaled range 
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function is defined by 

i?( fe ) /s maxi< m<n (y^ fc) (n)) -mini< m<n (y^ fc) (n)) 

-77rn) = —. v ' - ~ =± '- k = l,...,l. (2) 



1 For completeness let us point out that early instrumental or paleoclimatic data can be used 
to extend the North Atlantic Oscillation index back to 1823 or even 1675 25-28 . 



The average of the rescaled range in all boxes with an equal size n is next obtained 
and denoted by < R/S >. The above computation is then repeated for different values 
of n to provide a relationship between < R/S > and n, - which is expected to be a power 
law < R/S >~ n H if some scaling range and law exist; H is called the Hurst exponent. If 
H = 0.5, the signal is uncorrelated (white noise); the "walk" is like a Brownian motion. 
If H < 0.5, the signal is anticorrelated (blue noise); if H > 0.5, there are positive 
correlations in the signal (red noise). From Fig. 3 it is found that H = 0.55 ± 0.02 for 
the NAO index variations from 10 to 300 months (about 25 years). The departure from 
strict linearity on a log-log plot is usually attributed to too small box sizes, or to some 
periodic mode not finally taken into account through the assumed constant base line in 
Eq. (1); see discussions in references quoted here above and in the introduction. In order 
to test the robustness of the result we have checked the scaling properties of NAO index 
data series that are shorter than the original one by 5% (107 data points). In both cases 
when we delete the first 107 data points or the last 107 data points, the R/S analysis 
exponent has roughly the same value. The same goes true (not shown) for surrogate data 
series, i.e. when amplitudes are randomly displaced or multiplied by a random sign. It is 
certain that the H value near 0.55 indicates a weak deviation of the signal from Brownian 
motion. 

C. The detrended fluctuation analysis 

The detrended fluctuation analysis method has been recently much used in the mete- 
orological field 19 ' 30-32 . The method has the advantages of avoiding (seasonal-like) trends 
and non stationarity effects, intrinsic to the finite size of the data. The method consists 



in dividing the time series X = {xi,i = 1,...,N} into I boxes of equal size n. In the 
~k th box, the cumulative sum Y^ can be calculated as above, in the so called the first 
order DFA. Let Yf iti (n), be the best linear fit to the data in the k th box. The detrended 
fluctuation function is next calculated, i.e. dropping the (n), (f>\ — Y^' —Yr iti . The root 
mean square fluctuation is then given by 
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If the values of the time series are correlated, there is a power-law relationship between 
F{n) and n: F(n) ~ n a . A departure from linearity on a log-log plot, and the existence 
of crossovers (hereby one occurs near 240 months, see Fig.4) has been discussed 17 . Fig. 
4 shows that a = 0.54 ± 0.02 below 240 months, for NAO index variations. Again 
the robustness of the result is confirmed by analysing shorter or surrogate data series. 
Various considerations 15 indicate that a should be equal to H . Thus we can conclude 
that the value of the scaling exponent a is roughly the same as the one obtained within 
the rescaled range analysis. Notice that the DFA method, is clearly leading to extend 
the scaling properties of the NAO index toward smaller scales, less than 10 months as is 
found for the R/S analysis. 

The above findings confirm the existence of non trivial correlations (since a ^ 0.5, 
even though it is close to 0.50) within precise interval time ranges. They point out to the 
existence of physical phenomena described as fractional Brownian motions (Mandelbrot 
1982) thus with a fractal- like hierarchy of time scales. The result of a values larger than 
0.5 can be interpreted again through a persistence effect in the fluctuations 15 ' 29 . 



D. The Fokker-Planck equation 

In view of the above it is of interest to search whether these weak and persistent 
correlations can be found through the solution of a phenomenological evolution equation, 
like the Fokker-Planck equation 34 . Thus we focus on the variations Ax of the elements of 
the NAO series and the more so on their distribution in time. In order to do so we follow 
the method of Friedrich et al. 35 and reproduce it almost in extenso here below, for the 
technique is not necessarily familiar to most readers. 

In order to characterize the statistics of NAO changes, increments Axi, Ax 2 , ... for 
delay times Ati, At 2 , ••• at the same time t are considered. This leads to a set of 
p(Axi,Ati). Next the joint probability density functions are evaluated for various time 
delays Ati > At 2 > A^ > ... directly from the given data set, e.g. p(Axi, Ati, Ax 2 , At 2 ). 
Of course if two increments i.e. Ax\ and Aa; 2 are statistically independent, the joint pdf 
should factorize into a product of two probability density functions: 

p(Axi, Ati; Ax 2 , At 2 ) = p{Ax±, At 1 )p(Ax 2 , At 2 ). (4) 

leading to an isotropic single hill landscape in the Axi, Ax 2 plane. 

A complete characterization of the statistical properties of the data set in general 
requires the evaluation of joint pdVs p N (Axi,Ati;...;Ax^,At^) depending on N variables 
(for arbitrarily large N). In the case of a Markov process (a process without memory but 
governed by probabilistic conditions), an important simplification arises: The iV-point pdf 
p N is generated by the mere product of conditional probabilities p(Ax{ + i, Ati + i\Axi, Ati) 
itself equal to p(Axi + i, Ati + i; Axi, Ati)/p(Axi, Ati) for % = 1, ..., N — 1. The conditional 
probability is given by the probability of finding Axi + i values for fixed Ax,. As a necessary 



condition of Markov processes, the Chapman-Kolmogorov equation in its integral form 
reads 

p(Ax 2 , At 2 \A Xl , Ah) = f d(Ax i )p(Ax2,At 2 \Ax i ,At i )p(Ax i ,At i \Ax 1 ,At 1 ) (5) 

and should hold for any value of Ati, with Ai 2 < Ati < Ati, see Appendix A for a 
discussion in particular concerning large (and thus rare) events. As is well known, such 
a Chapman-Kolmogorov equation yields an evolution equation for the change of the con- 
ditional distribution functions p (Ax, At\Ax±, Ati) and p(Ax,At ) across the scales At. 
The Chapman-Kolmogorov equation when formulated in differential form yields a mas- 
ter equation, which can take the form of a Fokker-Planck equation 35-37 ' 24 . It is useful to 
use reduced time units, like r = log 2 (16/At), 2 

±p(Ax, r) = l-^D^(Ax, r) + ^^(Ax, r)]p(Ax, r) (6) 

in terms of a drift D^(Ax,r) and a diffusion coefficient D^ 2 \Ax,r) (thus values of r 
represent At,, % = 1, ....) Their functional dependence can be estimated directly from 
the moments M^ k > (known as Kramers- Moyal coefficients) of the conditional probability 
distributions: 

M {k) = -!- / dAx'(Ax - Ax) k p(Ax, r + Ar| Ax, r) (7) 

Ar J 

for different small Ar's, such that for Ar — ► 0, 

D (k \Ax,T)~him Ar ^ M (k) . (8) 



2 Why " 16" is chosen to be the normalizing value will be made clear below, but it has obviously 
not much effect at this stage. 
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After calculating such moments from the conditional probabilities, we find (Fig. 5) that 
the coefficient M^ shows a linear dependence for small Ax, while M^ can be approxi- 
mated by a polynomial of degree two in Ax. Therefore the type of fluctuation probability 
drift term D^ is well approximated by a linear function of Ax, whereas the diffusion 
term D^ 2 ' follows a function quadratic in Ax. For very large values of Ax the statistics 
becomes poorer and the uncertainty increases. 

From a careful analysis of the data based on the functional dependences of M ^ and 
7Vf( 2 ) (Fig. 5 (a-b)), the following approximations hold true: 

DM = -0.52 Ax + 0.04 for I Ax |< 5 (NAO units) 

(9) 

L> (2) = \ (0.21 Ax 2 - 0.02Ax + 4.2) for | Ax |< 5 (NAO units) 

Notice the range of validity of the simple analytical forms, thus the limit found for what 
would be called 38 "extreme events" or "outliers". Also observe that except for the in- 
dependent term, 2D^ ~ (D^) 2 , the strict equal sign being a request for indicating an 
absolute lack of intermittency in turbulence 39-42 . 

The FPE for the distribution function is known to be equivalent to a Langevin equation 
for the variable, i.e. Ax here, within the Ito interpretation 24 ' 34 

d 



—Ax(r) = D (1) (Ax(r),r)+n(r)JD( 2 )(Ax(r),r), (10) 

dr 

where n(r) is a fluctuating 5-correlated force with Gaussian statistics, i.e. < n(r) n(r')> = 
2S(t—t'). An interpretation of the analogy between these drift and diffusion coefficients 
and those usually employed in fluid mechanics is given in Appendix. It may be worthwhile 
to emphasize here that (i) a negative slope value for D^> indicates a sort of restoring or 
damping force for the evolution of Ax; (ii) the observed quadratic dependence of the 
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diffusion term D^ is essential for the logarithmic scaling of the intermittence parameter 
in turbulence 41,42 . 

Using those analytical expressions of the empirically derived Kramers-Moyal coeffi- 
cients, Eq.(6) can now be integrated, thereby leading to a test of the Markov process 
assumption. For the above change of variables, r = log2(lQ/At), we take the observed 
distribution of the time series for At = 16 months thus at "large time", as the initial 
condition for the integration, - considering that there is no propagation of anomalous 
correlation to be expected after such a time lag. Indeed the scaling ranges of H and a 
found here above (Figs. 3 and 4) indicate that the correlations are quasi identical for 
At = 16 months or At = 240 — 400 months, whence justifying the normalization "16" 
as the safest lowest boundary range for quasi Brownian fluctuations. Starting that far 
with a Gaussian as initial condition, the results of the integration (for r = 0,1,2,3,4, or 
At= 16, 8, 4, 2, 1) are rather trivial. The variance and the mean are found to vary very 
weakly. The experimental data pdf's are shown in Fig. 6, - together even with the pdf 
for At = 400 months, - the latter being indistinguishable from the At = 16 months case. 
It is readily remarkable that the pdf's are close to the Gaussian form in the interesting 
NAO (Ax) index range, - with some marginal deviation for the extreme (and rare) events. 

E. Discussion 

It is therefore confirmed that the NAO is a complex phenomenon which is almost 
Markovian. This stresses the need to insert appropriate feedback mechanisms into any 
model evolution equation(s), with an appropriate (red) noise term. This has been recently 
discussed 4,5 ' 43 ' 20 . Long-range fractionally integrated noise seems indeed to provide a better 
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fit of the NAO SLP wintertime index over the period 1864 - 1998 than does either 
stationary red noise or a non- stationary random walk 9 . 

The persistence of the NAO index fluctuations, i.e. SLP fluctuations, is in agreement 
with the persistence of the sea surface temperature fluctuations at different sites in the 
North Atlantic as found by Monetti et al. 44 . Some reasons for the above can be found in 
studies based on circulation-like models 45 . 

III. CONCLUSIONS 

In summary, the aims of this paper are twofold : (i) to search whether scaling ranges 
exist in the North Atlantic Oscillation pressure index ; (ii) to derive its FPE and check 
the validity of the Markovianity assumption. This allows one to examine different time 
scales on the same footing, - a fundamental need in geophysics 46 . 

It is found that the lack of departure from a Gaussian process definitely is a new 
quality of the NAO index data set, - not perceptible with a rescaled range analysis, or a 
DFA, as done above, nor with spectral studies. 

However, it seems that the actual NAO index is not very useful 3 . Thus one might have 
to request other measurements for better predictability of climate and weather in Europe 
and the Northern Hemisphere, e.g. at other locations. We might also suggest studies on 
"not-monthly-averaged" indices, - a monthly average being strangely unphysical in our 



3 This is in agreement with conclusions from a recent paper casting doubt on the NAO-global 
warming relationship 47 , and indicating a too strong influence of the Azores data with respect 
to the Iceland one. See also Czaja and Frankignoul 48 
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opinion. 

Yet it is emphasized that the FPE provides the complete knowledge as to how the 
statistics of correlations in the index distribution change on different delay times. Since 
this includes an analysis in time t for a scalar Ax, it seems that the findings could 
be implemented in atmospheric weather low dimensional vector — models 49,50 . Further 
work in line with the above should be to relate the FPE to an analytical solution, e.g. 
with a model of the turbulence-like dynamics as was done for financial indices 51 or ionic 
transport through membranes 52 through a Beck-Tsallis 53 ' 54 nonextensive thermodynamics 
approach. 
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Appendix 

The coefficients called drift and diffusion used in the main text pertain to the evolu- 
tion of the pdf; they are usually describing the motion of particles in fluid mechanics. For 
example the diffusion coefficient occurs in the Einstein or Langevin equation of Brownian 
motion, as 

D = ^L (if) 

6nr]a 

where T is the bath temperature, n the dynamic viscosity of the fluid and a the diameter 
of the particle, such that the evolution of the particle is described as a function of time t 
by < x 2 > = 2Dt, solution of the Langevin equation 

M d ^ = -6n V a^ + R(t), (12) 

where M is the mass of the particle and R(t) a random force with zero mean. 

The Langevin equation is equivalent to the standard diffusion equation for a proba- 
bility density 

dpjx, t) ^ d 2 p{x,t) 

for which the solution is a Gaussian 

p(x,t) = -=± e~ x2 ^ Dt \ (14) 

The FPE written in the main text contains an extra term to Eq.(14); let it be rewritten 
here as 

j t p(x,t) = [- J^«(x,t) + -JL D W(x,t)]p(x,t), (15) 
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from which the mean < x > and the variance < a 2 > can be defined as usual. Keeping 
only the linear term in D^\ with a coefficient D[ and the independent Dq and quadratic 
D<2 terms in D^ 2 \ from Eq.(9), one easily obtains the evolution of the "particle" as 

(x(t)) = (x (t))e^ )+D ^ (16) 

(2) 

and a similar equation for the variance, but also containing D , from which one observes 
that D^ and D^ are true drift and diffusion coefficients. Note the time scale given by 
the inverse of D^ and D^ 2 \ i.e., about 1 month. 

Knowing that Ax is the NAO index, a difference in pressure, (AP) we can roughly 
rewrite the "official" diffusion coefficient, Eq.(12), as 

< (AP) 2 >= f^t, (17) 

and "interpret it", suggesting that in further and more precise work, one could develop 
a model relating the changes in pressure (between Iceland and the Azores) with a tem- 
perature field (in principle a temperature gradient, rather than the mean temperature of 
the bath). 

No need to say that the solution of a Brownian particle in a (rotating) bath under a 
temperature gradient and with a noise force term is indeed what a good weather model 
is (or should be). 



16 



REFERENCES 

1 Bhattacharya, J., and P.P. Kanjilal, 2000: Revisiting the role of correlation coefficient 
to distinguish chaos from noise, Eur. Phys. J. B, 13, 399-403. 

2 Provenzale, A., L.A. Smith, R. Vio, and G. Murante, 1992: Distinguishing between 
low-dimensional dynamics and randomness in measured time series, Physica D, 58, 
31-49. 

3 Rowlands, G., and J.C. Sprott, 1992: Extraction of dynamical equations from chaotic 
data, Physica D, 58, 251-259. 

4 Penland, C, and L. Matrosova, 1998: Prediction of tropical Atlantic sea surface tem- 
peratures using linear inverse modeling, J. Climate, 11, 483-496. 

5 Penland, C, and P.D. Sardeshmukh, 1995: The optimal growth of tropical sea surface 
temperature anomalies, J. Climate, 8, 1999-2024. 

6 Marshall, J., Y. Kushnir, D. Battisti, P. Chang, A. Czaja, R. Dickson, M. McCartney, 
R. Saravanan, and M. Visbeck, 2001: North Atlantic Climate Variability: phenomena, 
impacts and mechanisms, Inter. Jour. Climat., 21 1863-1898 

7 Hurrell, J.W., Y. Kushnir, and M. Visbeck, 2001: The North Atlantic Oscillation, 
Science, 291, 603-604. 

8 Palmer, T.N. , 2000: Predicting uncertainty in forecasts of weather and climate, Rep. 
Prog. Phys., 63 71-116. 

9 Greatbatch, R. J., 2000: The North Atlantic Oscillation, Stochastic Environmental 
Research and Risk Assessment, 14, 213-242. 

17 



10 Wanner, H., S. Bronnimann, C. Casty, D. Gyalistras, J. Luterbacher, C. Schmutz, D.B. 
Stephenson, and E. Xoplaki : 2001, North Atlantic Oscillation - concepts and studies, 
Surv Geophys , 22, 321-382. 

11 Wallace, J.M., 2000: North Atlantic Oscillation / Annular Mode: Two paradigms ? 
One Phenomenon, Quart. J. Royal Met. Soc, 126, 791-805. 

12 Ambaum, M. H. P., B. J. Hoskins, and D. B. Stephenson, 2001: Arctic Oscillation or 
North Atlantic Oscillation?, J. Climate, 14, 3495-3507. 

13 Stephenson, D.B. , V. Pavan, and R. Bojariu, 2000: Is the North Atlantic Oscillation 
a random walk?, Int. J. Climatoi, 20, 1-18. 

14 Hurst, H.E., 1951: Long-term storage capacity of reservoirs, Trans. Am. Soc. Civil 
Eng., 116, 770-808. 

15 Turcotte, D.L., 1999: Fractal and Chaos in Geology and Geophysics, Cambridge Univ. 
Press, Cambridge, UK. 

16 Azbel, M.Ya., 1995: Universality in a DNA Statistical Structure, Phys. Rev. Lett., 75, 
168-171. 

17 Hu, K., P. Ch. Ivanov, Zhi Chen, P. Carpena, and H. E. Stanley , 2001: Effect of trends 
on detrended fluctuation analysis, Phys. Rev. E, 64, 11 114 - 11 121. 

18 Pelletier, J.D., 1997: Analysis and modeling of the natural variability of climate, J. 
Climate, 10, 1331-1342. 

19 Talkner, K., and R.O. Weber, 2000: Power spectrum and detrended fluctuation analy- 
sis: Application to daily temperatures, Phys. Rev. E , 62, 150-160. 

18 



20 Fernandez, I., C.N. Hernandez, and J.M. Pacheco, 2003: Is the North Atlantic Oscilla- 
tion just a pink noise?, Physica A, 323, 705-714. 

21 Hasselmann, K., 1976: Stochastic climate models. Part 1: Theory, Tellus, 28, 473-485 

22 Frankignoul, C, and K. Hasselmann, 1977: Stochastic climate models. Part II: ap- 
plication to sea-surface temperature anomalies and thermocline variability, Tellus, 29, 
289-305. 

23 Leith, C.E., 1978: Predictability of climate, Nature, 276, 352-355. 

24 Risken, H., 1984: The Fokker-Planck Equation, Springer- Verlag, Berlin. 

25 Jones, P.D., T. Jonsson, and D. Wheeler, 1997: Extension to the North Atlantic Oscil- 
lation using early instrumental pressure observations from Gibraltar and South- West 
Iceland, Int. J. ClimatoL, 17, 1433-145. 

26 Jones, P.D., K.R. Briffa, T.P Barnett, and S.F.B. Tett, 1998: High resolution paleocli- 
matic records for the last millenium: interpretation, integration and comparison with 
general circulation model control- run temperatures, The Holocene, 8, 455-471. 

27 Luterbacher, J., C. Schmutz, D. Gyalistras, E. Xoplaki, and H. Wanner, 1999 : Recon- 
struction of monthly NAO and EU indices back to AD 1675, Geophys. Res. Lett., 26, 
2745-2748. 

28 Schmutz, C, J. Luterbacher, D. Gyalistras, E., Xoplaki, and H. Wanner, 2000 : Can 
we trust proxy-based NAO index reconstructions?, Geophys. Res. Lett., 27, 1135-1138. 

29 Malamud, B.D., and D.L. Turcotte, 1999: Self-amne time series: Measures of weak and 
strong persistence, J. Stat. Plann. Infer. 80, 173-196. 

19 



30 Koscielny-Bunde, E., A. Bunde, S. Havlin, H. E. Roman, Y. Goldreich, and H.-J. 
Schellnhuber, 1998: Indication of a universal persistence law governing atmospheric 
variability Phys. Rev. Lett, 81, 729-732. 

31 Ivanova, K., and M. Ausloos, 1999: Application of the Detrended Fluctuation Analysis 
(DFA) method for describing cloud breaking Physica A , 274, 349-354. 

32 Ivanova, K., M. Ausloos, E.E. Clothiaux, and T.P. Ackerman, 2000: Break-up of stra- 
tus cloud structure predicted from non-Brownian motion liquid water and brightness 
temperature fluctuations, Europhys. Lett., 52, 40-46. 

33 Mandelbrot, B.B., 1982: The Fractal Geometry of Nature, Freeman, San Francisco. 

34 Pecseli, H.L., 2000: Fluctuations in Physical Systems, Cambridge Univ. Pres, Cam- 
bridge, UK. 

35 Friedrich, R., J. Peinke, and Ch. Renner, 2000: How to quantify deterministic and 
random influences on the statistics of the foreign exchange market, Phys. Rev. Lett., 
84, 5224-5228. 

36 Hanggi, P., and H. Thomas, 1982: Stochastic processes : time evolution, symmetries 
and linear response, Phys. Rep., 88, 207-319. 

37 Gardiner, C.W., 1983: Handbook of Stochastic Methods, Springer- Verlag, Berlin. 

38 Grasso, J.R., and D. Sornette, 1998: Testing self-organized criticality by induced seis- 
micity, J. Geophys.Res., 103, 29965-29987. 

39 Kolmogorov, A.N., 1941a: Local structure of turbulence in an incompressible viscous 
fluid for very large Reynolds numbers, Doklady Akad. Nauk SSSR, 30, 299-303. 

20 



40 Kolmogorov, A.N., 1941b: On degeneration of isotropic turbulence in an incompressible 
viscous liquid, Doklady Akad. Nauk SSSR, 31, 538-541. 

41 Parisi, G., and U. Frisch, 1985: in Turbulence and Predictability in Geophysical Fluid 
Dynamics and Climate Dynamics, Ghil M., Benzi R. and Parisi G., Eds. North Holland, 
New York. 

42 Frisch, U., 1995: Turbulence, Cambridge Univ. Press, Cambridge UK. 

43 Johnson, S.D., D.S. Battisti, and E.S. Sarachik, 2000: Empirically Derived Markov 
Models and Prediction of Tropical Pacific Sea Surface Temperature Anomalies, J. Cli- 
mate, 13, 3-17. 

44 Monetti, R.A., S. Havlin, and A. Bunde, 2003: Long-term persistence in the sea surface 
temperature fluctuations, Physica A, 320 581-589. 

45 Josey, S. A., E.C. Kent, and B. Sinha, 2001 : Can a State of the Art Atmospheric 
General Circulation Model Reproduce Recent NAO Related Variability at the Air-Sea 
Interface?, Geophys. Res. Lett, 28, 4543-4546. 

46 Fraedrich, K. and C.-D. Schonwiese, 2002: Space-Time Variability of the European 
Climate in The Science of Disasters, A. Bunde, J. Kropp, H.J. Schellnhuber (eds.), 
Springer, Berlin, pp. 105-139. 

47 Jonsson, T. and M. Miles, M. 2001: Anomalies in the seasonal cycle of sea level pressure 
in Iceland and the North Atlantic Oscillation, Geophys. Res. Lett., 28, 4231-4234. 

48 Czaja, A., and C. Frankignoul, 2002: Observed impact of Atlantic SST anomalies on 
the North Atlantic Oscillation, J. Climate, 15, 606-623. 



21 



49 Grabowski, W.W., and P.K. Smolarkiewicz, 1999: CRCP: a cloud resolving convection 
parameterization for modeling the tropical convecting atmosphere, Physica D, 133, 
171-178. 

50 Ragwitz, M., and H. Kantz, 2000: Detecting non-linear structure and predicting tur- 
bulent gusts in surface wind velocities, Europhys. Lett., 51, 595-601. 

51 Ausloos, M. and K. Ivanova, 2003: Dynamical model and nonextensive statistical me- 
chanics of a market index on large time windows, Phys. Rev. E , 68, 046122. 

52 Ausloos, M., K. Ivanova and Z. Siwy, 2004: Searching for self-similarity in switch- 
ing time and turbulent cascades in ion transport through a biochannel. A time delay 
asymmetry, Physica A , 336, 319-333. 

53 Tsallis, C, 1988: Possible generalization of Boltzmann-Gibbs statistics J. Stat. Phys., 
52, 479-. 

54 Beck, C, 2001: Dynamical foundations of nonextensive statistical mechanics, Phys. 
Rev. Lett. 87, 180601. 



22 



Figure Captions 

Figure 1 Time evolution of the monthly averaged NAO index fluctua- 
tions from January 1825 to November 2002 (2135 points) available on the 
web sites http : / / www.cru.uea. ac.uk/ ftpdata/nao.dat and updated at http : 
// www.cru.uea.ac.uk/limo/projpages/nao-update.htm) 

Figure 2 Integration of the NAO index fluctuations signal shown on Figure 1 

Figure 3 Rescaled range analysis of the NAO index signal. Inserts: Rescaled range 
analysis of the NAO index signal shortened by 5% at the beginning of the data series (left 
upper panel) or at the end (right lower panel) 

Figure 4 Detrended fluctuation function of the (integrated) NAO index signal. In- 
serts: Detrended fluctuation function of the (integrated) NAO index signal after shorten- 
ing by 5% at the beginning of the data series (left upper panel) or at the end (right lower 
panel) 

Figure 5 Kramers-Moyal coefficients (a) M^ and (b) M^ estimated from the 
empirical conditional density probability of the NAO distribution values. The solid curves 
represent a linear and a quadratic fit, respectively, excluding large events 

Figure 6 Raw data (symbols), and theoretical (Gaussian) pdf (solid line) comparing 
the NAO fluctuation distribution functions for various time lags, At = 16, 8, 4, 2, 1 
months (or for r = 0,1,2,3,4). The case At = 400 months is also "shown" for comparison 
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